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Analyzing power in elastic scattering of 6 He from polarized proton target at 
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The vector analyzing power has been measured for the elastic scattering of neutron-rich 6 He from 
polarized protons at 71 MeV/nucleon making use of a newly constructed solid polarized proton 
target operated in a low magnetic field and at high temperature. Two approaches based on local 
one-body potentials were applied to investigate the spin-orbit interaction between a proton and a 
6 He nucleus. An optical model analysis revealed that the spin-orbit potential for 6 He is characterized 
by a shallow and long-ranged shape compared with the global systematics of stable nuclei. A semi- 
microscopic analysis with a a+n+n cluster folding model suggests that the interaction between 
a proton and the a core is essentially important in describing the p+ 6 He elastic scattering. The 
data are also compared with fully microscopic analyses using non-local optical potentials based on 
nucleon-nucleon g-matrices. 

PACS numbers: 24.10.Ht, 24.70.+S, 25.40.Cm, 25.60.Bx, 29.25.Pj 



I. INTRODUCTION 

Spin-orbit coupling in atomic nuclei is an essential fea- 
ture in understanding any reaction and nuclear struc- 
ture related to it. One of the direct manifestations of 
that spin-orbit coupling in nuclear reactions, is the po- 
larization phenomenon in nucleon elastic scattering [1- 
3], Characteristics of the spin-orbit coupling between 
a nucleon and stable nuclei have been well established 
by analyses of measured vector analyzing powers in the 
elastic scattering of polarized nucleons on various targets 
over a wide range of incident energies [4-7] . 

On the other hand, the spin-orbit coupling of a nu- 
cleon with unstable nuclei might be considerably differ- 
ent from that with the stable nuclei. Some neutron-rich 
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nuclei with small binding energies arc known to have very 
extended neutron distributions [8]. Since the spin-orbit 
coupling is essentially a surface effect, it is natural to ex- 
pect that the diffused density distribution of a neutron- 
rich nucleus may significantly effect the radial shape and 
depth of the spin-orbit potential. The purpose of this 
work is to investigate the characteristics of the spin-orbit 
potential between a proton and 6 He; a typical neutron- 
rich nucleus. 

Experimental determination of the spin-orbit poten- 
tial strongly owes to measurements and analyses of the 
vector analyzing powers. However, until recently, ana- 
lyzing power data were not obtained in the scattering 
which involves unstable nuclei. This was mainly due to 
the lack of a polarized proton target that is applicable to 
radioactive ion (RI) beam experiments. Rl-beam exper- 
iments induced by light ions are usually carried out un- 
der inverse-kinematics conditions, where energies of recoil 
protons can be as low as 10 MeV. Conventional polarized 
proton targets [9, 10], based on the dynamic nuclear po- 
larization method, require a high magnetic field and low 
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temperature such as 2.5 T and 0.5 K, respectively. It is 
impossible to detect the low-energy recoil protons with 
sufficient angular resolution under these extreme condi- 
tions. For the application in Rl-beam experiments, we 
have constructed a solid polarized proton target which 
can be operated under low magnetic field of 0.1 T and 
at high temperature of 100 K [11-15]. The electron po- 
larization in photo-excited aromatic molecules is used to 
polarize the protons [16, 17]. A high proton polarization 
of about 20% can be achieved in relatively "relaxed" op- 
erating conditions described above, since the magnitude 
of the electron polarization is almost independent of the 
magnetic field strength and temperature. 

We have measured the vector analyzing power for the 
p+ 6 He elastic scattering at 71 MeV/nucleon [18] using 
the solid polarized proton target, newly constructed for 
Rl-beam experiments. 6 He is suitable for the present 
study since it has a spatially extended distribution due 
to a small binding energy. In addition, from an experi- 
mental viewpoint, the p+ 6 He elastic scattering measure- 
ment is relatively easy to perform since 6 He does not 
have a bound excited state. This allows us to identify the 
elastic-scattering event only by detecting 6 He and a pro- 
ton in coincidence. The analyzing powers thus measured 
are the first data set that can be used for quantitative 
evaluation of the spin-orbit interaction between a proton 
and an unstable 6 Hc nucleus. The essence of these mea- 
surements has been published in Rcf. [18] together with 
two kinds of theoretical analyses by folding models; one 
assumes a fully antisymmetrized large-basis shell model 
for 6 He with the ^-matrix interaction and the other an 
a+n+n cluster model for 6 He with a p-n effective in- 
teraction and a realistic p-a static potential. The main 
purpose of the present paper is to give more details of 
the experiment and present an additional analysis of the 
experimental data using a one-body p- 6 He optical poten- 
tial. The analysis exhibits remarkable characteristics for 
the spin-orbit part of that potential. Then it becomes 
important to investigate if such a potential can be de- 
rived theoretically from any model of 6 Hc. As the first 
approach we examined the a+n+n folding potential in 
more detail, since important contributions of the a clus- 
ter are suggested by the fact that the measured A y for 
6 He is similar to that for 4 He [18], when plotted versus the 
momentum transfer of the scattering. To identify effects 
of the clusterization, we also calculated the p- 6 Hc folding 
potential for a 2p+An non-cluster model of 6 He and com- 
pared the results with those of the a+n+n cluster model. 
Hereafter, they are referred to as the ann cluster folding 
(CF) model and nucleon folding (NF) one, respectively. 
In addition, the data are also compared with fully mi- 
croscopic calculations using non-local optical potentials. 
In this model, non-locality of the p- 6 He interaction, a 
consequence of the Pauli principle leading to nucleon ex- 
change scattering amplitudes, is taken into account ex- 
plicitly. Three sets of single-particle wave functions, as 
well as the required one-body density matrix elements 
determined from a large-basis shell model for 6 He, have 



been used in these calculations. 

The present paper is subdivided as follows. In Section 
2, details of the experimental method are described. In 
Section 3, the method of the data reduction is presented. 
Section 4 deals with the phenomcnological optical model 
analysis. Section 5 is devoted to the details of the ann 
cluster folding calculation and the nucleon folding cal- 
culation. In Section 6, the data are compared with the 
analysis by the non-local g-folding optical potentials. Fi- 
nally, a short summary of the obtained results is given in 
Section 7. 



II. EXPERIMENT 
A. Experimental setup 

The experiment was carried out at the RIKEN Ac- 
celerator Research Facility (RARF). The 6 He beam was 
produced through the projectile fragmentation of a 12 C 
beam with an energy of 92 MeV/nucleon bombarding a 
primary target. As that primary, we used a rotating 9 Be 
target [19] to avoid heat damage by the beam. A thick- 
ness of the target was 1480 mg/cm 2 . The 6 He particles 
were separated by the RIKEN Projectile-fragment Sepa- 
rator (RIPS) [20] based on the magnetic rigidity and the 
energy loss of fragments. The energy of the 6 He beam 
was 70.6±1.4 MeV/nucleon at the center of the secondary 
target. The purity of the beam was 95%. 

The solid polarized proton target was placed at the 
final focal plane of RIPS. Figure 1 illustrates the ex- 
perimental setup of the target and detectors. The most 
prominent advantage of the target is its relaxed operat- 
ing conditions, i.e. a low magnetic field of 0.1 T and high 
temperature of 100 K. These conditions allow us to de- 
tect recoil protons whose energies are as low as 10 MeV. 
Details of the target will be described in the following 
subsection. 

A detector system consisted of two subsystems: one 
for scattered particles and the other for recoil protons. 
Detection of the recoil protons with energies as low 
as 10 MeV is essential for the selection of the clastic- 
scattering events. The scattering angles of protons were 
determined by single- wire drift chambers (SWDC). The 
SWDCs were placed 138.5 mm away from the target on 
both left and right sides of the beam axis as shown in 
Fig. 1. They covered an angular region of 39°-71° (hor- 
izontal) and ±9.7° (vertical) in the laboratory system. 
Their position resolution and detection efficiency were 
found to be 2.6 mm (FWHM) and 99.3%. For the mea- 
surement of the total energy of protons, we used CsI(Tl) 
scintillation detectors. They were placed just behind 
the SWDCs. Light output from the CsI(Tl) crystal was 
detected by photo-multiplier tubes. The front side of 
the CsI(Tl) scintillator was covered by the thin carbon- 
aramid film with a thickness of 12 ^m. Material thickness 
of the film, the SWDC, and air between the detectors was 
24 mg/cm 2 in total. Energy loss of 10 MeV protons in 
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FIG. 1: (Color online) Experimental setup of the secondary target and detectors is shown. 



these materials is 1.2 MeV, which does not prevent the 
detection. 

A multi-wire drift chamber (MWDC) was used to re- 
construct the trajectories of scattered particles. Scatter- 
ing position on the secondary target was determined from 
the reconstructed trajectory. The MWDC was placed at 
880 mm downstream of the target. It has a sensitive area 
of 640 mm (horizontal) x 160 mm (vertical) and covered 
an angular region of ±16° x ±4° in the laboratory system. 
The configuration of the planes of the MWDC is X-Y- 
X'-Y'-X'-Y'-X-Y, where "X(Y)-plane" has anode-wires 
oriented along the vertical (horizontal) axis. The planes 
with primes are displaced with respect to the "unprimed" 
planes by half the cell size. The cell size is 20 mmx20 mm 
for the X-plane and 10 mmx 10 mm for the Y-plane. The 
material of the anode wire is gold-plated tungsten with 
a diameter of 30 fim. Negative high voltages were ap- 
plied to the cathode and potential wires: —2.85 kV for 
the X (X')-planes and -2.15 kV for the Y (Y')-planes. A 
gas mixture of Ar (50%) and C2H 6 (50%) was used. Po- 
sition resolution and detection efficiency of the MWDC 
were found to be 0.2 mm (FWHM) and 99.8%. For iden- 
tification of scattered particles, we used a plastic scin- 
tillation detector array placed just behind the MWDC. 
The first and second layers with thicknesses of 5 mm and 
100 mm provided information of the energy loss and the 
total energy of scattered particles. The total number of 
the beam particles was counted with a beam monitor 
placed between the secondary target and the MWDC. A 
50 mm ff x 50 mm^ x 10 mm D plastic scintillator was used 
for the beam monitor. A beam stopper made of a copper 



block was placed just behind the beam monitor. 



B. Solid polarized proton target 

The solid polarized proton target, used in the measure- 
ment, can be operated in a low magnetic field of 0.1 T and 
at a high temperature of 100 K. These relaxed operation 
conditions allow us to detect low-energy recoil protons 
without losing angular resolution. This capability is in- 
dispensable to apply the target to scattering experiments 
carried out under the inverse kinematics condition. The 
proton polarization of about 20% has been achieved [15] 
under such relaxed conditions by introducing a new po- 
larizing method using electron polarization in triplet 
states of photo-excited aromatic molecules [16, 17]. A 
single crystal of naphthalene (CioHg) doped with a small 
amount of pentacene (C22H14) is used as the target ma- 
terial. Protons in the crystal are polarized by repeating 
a two-step process: production of electron polarization 
and polarization transfer. In the first step, pentacene 
molecules are optically excited to higher singlet states. 
A small fraction of them decays to the first triplet state 
via the first excited singlet state by the so-called inter- 
system crossing. Here, electron population difference is 
spontaneously produced among Zecman sublevels of the 
triplet state [16]. In the second step, the electron popu- 
lation difference between two Zeeman sublevels, namely 
electron polarization, is transferred to the proton polar- 
ization by the cross-relaxation technique [17]. 

As the target material, we used a single crystal of naph- 
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thalenc doped with 0.005 mol% pentacenc molecules. 
The crystal was shaped into a thin disk whose diame- 
ter and thickness are 14 mm and 1 mm (116 mg/cm 2 ), 
respectively. The number of hydrogens per unit area was 
4.29±0.13xl0 21 /cm 2 . In order to reduce the relaxation 
rate, the target crystal was cooled down to 100 K in 
a cooling chamber with the flow of cold nitrogen gas. 
The cooling chamber was installed in another chamber 
as shown in Fig. 1. Heat influx to the cooling chamber 
was reduced by the vacuum kept in the intervening space 
between these two chambers. Each chamber has one win- 
dow (6 /im-thick Havar foil) on the upstream side for the 
incoming Rl-beam, two glass windows for the laser irra- 
diation, and three windows (20 ^m-thick Kapton foil) on 
the left, right, and downstream sides for the detection of 
recoil and scattered particles. 

A static magnetic field was applied on the target crys- 
tal by a C-type electromagnet to define the polarizing 
axis. The gap and the diameter of the poles were 100 mm 
and 220 mm, respectively. The strength of the magnetic 
field in the present experiment was 91 mT; a value much 
higher than that of the crystal field (« 2 mT). While the 
effects of the magnetic field on the scattering angles of 
6 He particles and protons were sufficiently small (about 
0.07° and 0.2-0.8°, respectively), they were properly cor- 
rected in the data analysis. 

The target crystal was irradiated by the light of two 
Ar-ion lasers with a power of 25 W each in the multi- 
line mode. Wavelengths of main components of the light 
were 514.5 nm (10 W) and 488.5 nm (8 W). The laser 
light was pulsed by a rotating optical chopper. Typi- 
cally the pulse width and repetition rate were 12-14 jj,s 
and 1 kHz. Microwave (MW) irradiation and a magnetic 
field sweep are required in the cross-relaxation method. 
For the MW irradiation, the target crystal was installed 
in a resonator. In order to detect low-energy recoil pro- 
tons, we employed a thin cylindrical loop-gap resonator 
(LGR [21]) made of 25 /um-thick Teflon film. Copper 
stripes with a thickness of 4.4 ^m were printed on both 
sides of the film. The MW frequency was 3.40 GHz. The 
LGR was surrounded by a cylindrical MW shield made 
of 12 jum-thick aluminum foil. For the cross-relaxation, 
the magnetic field was swept from 88 mT to 94 mT at 
the rate of 0.36 mT/^s, simultaneously with the MW ir- 
radiation, by applying a current to a small coil placed in 
the vicinity of the target material. 

Proton polarization was monitored during the experi- 
ment by the pulse NMR method. A radio- frequency (RF) 
pulse with a frequency and a duration of 3.99 MHz and 
2.2 /is was applied to a 19 mm</> NMR coil covering the 
target crystal. The free induction decay (FID) signal was 
detected by the same coil. We carried out the absolute 
calibration to relate the FID signal to the proton polar- 
ization by measuring the spin-asymmetry in the p+ He 
elastic scattering. Details of the calibration procedure 
are described in Appendix. A. 

Devices located near to the target, namely the LGR, 
the MW shield, the field sweeping coil, and the NMR 



coil, were fabricated with hydrogen-free materials to pre- 
vent production of background events. Table I shows 
the material thicknesses of the devices that recoil pro- 
tons penetrate. Energy losses of the 20 MeV protons in 
these materials are sufficiently small for the detection as 
summarized in Tabic I. 

The target polarization during the experiment is shown 
in Fig. 2 as a function of time. The polarization was 
built up for the first 40 hours and reached the maxi- 
mum value of 20.4±3.9%. The target was then irradi- 
ated by a 71 MeV/nucleon 6 He beam for 55 hours, by 
a 80 MeV/nucleon 4 He beam for the following 25 hours, 
and again by the 6 He beam for 60 hours. The magni- 
tude of average polarization was found to be 13.8±2.7%. 
The target polarization slowly decreased as a function 
of time, which is due to beam-irradiation damage in the 
target material. This radiation damage increased the re- 
laxation rate of the target material from T = 0.127(6) h -1 
before the experiment to T = 0.295(4) h _1 after the beam 
irradiation. The direction of the target polarization was 
reversed three times during the measurement to cancel 
spurious asymmetries. The 180° pulse NMR method was 
used here. Reversal efficiency of 60-70% was achieved. 
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FIG. 2: Target polarization is shown as a function of time. 
The open squares, closed circles, and open circles indi- 
cate the target polarization during the polarization buildup, 
the p+ 6 He elastic-scattering measurement, and the p+ 4 He 
elastic-scattering measurement, respectively. 



III. DATA REDUCTION 

A. Data analysis 

In principle, elastic-scattering events of the 6 He from 
protons can be identified by the coincidence detection of 
scattered 6 Hc particles and recoil protons, since the 6 He 
does not have a bound excited state. Note that the first 
excited state of 6 He, which is the 2 + state at 1.87 MeV, is 
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Material 


Thickness (mg/cm^) 


Energy loss (MeV) 


Target crystal (Naphthalene) 


- 336 


- 9.5 


-LvjrK (ienon, Lu roil J 


n q 

y.o 


no n a 
\).Z — U.4 


Microwave shield (Al foil) 


3.2 


0.05 - 0.1 


Cooling gas (N2) 


13.5 


0.3 - 0.6 


Window (Kapton film) 


20 


0.5 - 1.0 


Total 


46 - 382 


1.1 - 11.6 



TABLE I: Thicknesses of the materials of target devices and energy losses of 20 MeV recoil protons in them. 



above the two-neutron breakup threshold (0.975 MeV). 
Thus, any excited 6 He particles decay into a+n+n sys- 
tems before reaching the detectors. 

Scattered particles were identified by the standard AE- 
E method. Figure 3 shows a two-dimensional plot of 
the total energies of scattered particles E versus their 
energy losses AE, where loci of tritons, 4 He, 6 He, and 8 Li 
are found. Tritons and 8 Li are the contamination in the 
secondary beam. Most of 4 He particles were produced by 
6 Hc dissociation in the secondary target. However, some 
originated from 6 Hc reactions in the plastic scintillators. 
So to count all of the p+ 6 He elastic-scattering events, the 
particle identification gate includes most of the 4 He locus 
as shown by solid curves in Fig. 3. The contribution of 
the dissociation reaction, which is not excluded by this 
gate, was subtracted using a kinematics relation. This is 
described after the response of the recoil proton detectors 
is considered. 
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FIG. 3: Two-dimensional plot of the total energies of scat- 
tered particles versus their energy losses. Solid curves indicate 
the particle-identification gate. 



Figure 4 shows a two-dimensional scatter plot of the 
proton energies versus their scattering angles in the 
center-of-mass system, 6 c . m .. The kinematic locus of the 
clastic scattering is clearly identified, while backgrounds 



from other reaction channels such as p( 6 He, p 4 He) are 
also evident. The kinematic locus of elastic-scattering 
events shows that the recoil protons were properly de- 
tected outside of the target. It should be noted that this 
correlation was not used for the event selection, since it 
would cause a loss of the events at forward angles. 
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FIG. 4: Two-dimensional plot of the proton energies versus 
their scattering angles. 



To discriminate elastic scattering from the background, 
we used the correlation of the azimuthal angles of pro- 
tons <j> p with those of scattered particles c/> sca tt.- In the 
case of the elastic scattering, a scattered 6 He and a recoil 
proton stay within a well defined reaction plane since the 
final state is a binary system. Thus, the difference of az- 
imuthal angles A<j) — <j) p — sca tt. makes a narrow peak at 
around 180°. This back-to-back correlation holds even if 
the scattered 6 He is dissociated in the plastic scintillator. 
In the case of other reactions, however, the azimuthal an- 
gle difference is more spread since their final states consist 
of more than two particles. 

Figure 5 shows the distribution of the azimuthal angle 
difference A0 fitted by a double- Gaussian function. The 
narrower component is reasonably identified as that of 
the elastic-scattering events. The peak width of 3.5° in 
sigma is consistent with the detector resolution of 3.1°. 
We selected the events of \A(f> - 180° | < 12.4°. The 
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background remaining in the gate was evaluated from 
the broader component and was subtracted. Contribu- 
tions of the inelastic scattering and other reactions such 
as breakup were removed in this way without losing the 
elastic-scattering yields. Figure 6 shows a background- 
subtracted two-dimensional plot of scattering angles in 
the center-of-mass system versus angles of scattered par- 
ticles. Center-of-mass scattering angles were deduced 
from recoil angles of the protons in the laboratory system, 
since the resolution of scattering angles of 6 He particles is 
insufficient due to the kinematic focusing. In Fig. 6, clear 
peaks of elastic-scattering events lie along the solid curves 
indicating the kinematics of the p+ 6 He elastic scattering. 
Small peaks at |#e Ho | re 4° originated from the ambigu- 
ity in the background subtraction. Yields of the p+ 6 He 
clastic scattering were obtained by counting the events 
of the elastic-scattering peaks in the typical width of 4° 
in 06 Ho- 
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FIG. 6: Scattering angle correlation between scattered par- 
ticles and recoil protons. Fhe solid curves indicate the kine- 
matics of the p+ 6 He elastic scattering. 



The present work demonstrates the applicability of the 
solid polarized proton target in the Rl-bcam experiment. 
The relaxed operation condition of the target, i.e. a low 
magnetic field of 0.1 T and high magnetic field of 100 K, 
enables us to detect the low-energy recoil protons. As 
described in the data analysis above, information on the 
trajectory of recoil proton is indispensable both in identi- 
fying the elastic-scattering events (Fig. 5) and in deduc- 
ing the scattering angle (Fig. 6). 



B. Experimental data 

The da/dVl of thep+ 6 He elastic scattering measured at 

71 MeV/nucleon are summarized in Table II. In the back- 
ward region, the uncertainty mainly results from statis- 
tics and from the ambiguity in the background subtrac- 
tion. In the forward angular region, 6 C m <C 60°, the main 
component of the uncertainty in da/dfl is the systematic 
uncertainty in the number of the incident particles (10%). 
The target was hit by only a fraction of beam particles 
since the size of the secondary beam was comparable to 
that of the target. The percentage of the beam particles 
incident on the target was determined from the beam 
profile and was found to be 65±7% of those counted by 
the beam monitor. The beam profile was measured with 
the MWDC by removing the beam stopper. Stability of 
the beam profile was confirmed by several measurements 
carried out before, during, and after the elastic-scattering 
measurement. 

The analyzing power A y is deduced with the standard 
procedure as 

where P denotes the target polarization. The values 
A r 's represent the yield of the clastic-scattering events 
where subscripts and superscripts denote the scatter- 
ing direction (left/right) and the polarization direction 
(up/down), respectively. The statistical uncertainty is 
expressed by 

AA y _ LR l~L 1 1 r~ 

This procedure allows us to minimize the systematic un- 
certainties originating from unbalanced detection efficien- 
cies and misalignment of detectors. The obtained A y are 
summarized in Table III. It must be noted that there is 
an additional scale error of 19% resulting from the un- 
certainty in the target polarization P (see Appendix. A). 

Figure 7 shows da/dQ. and A y for the p+ 6 He elastic 
scattering at 71 MeV/nucleon (closed circles: present 
work, open circles: Ref. [22]), those for the p+ 4 He at 

72 MeV/nucleon (open squares: Ref. [23]), and those for 
the p+ 6 Li at 72 MeV/nucleon (open triangles: Ref. [24]). 
The present data are consistent with the previous ones in 
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9cm. (deg 


) A0 c . m . (deg) 


% (mb/sr) 


AS ( mb /sr) 


42.1 


2.5 


5.02 


0.52 


47.1 


2.5 


2.03 


0.22 


52.1 


2.5 


0.796 


0.098 


57.4 


2.5 


0.454 


0.059 


62.3 


2.5 


0.360 


0.046 


67.3 


2.5 


0.226 


0.031 


72.3 


2.5 


0.172 


0.023 


77.3 


2.5 


0.127 


0.018 


82.2 


2.5 


0.064 


0.013 


87.2 


2.5 


0.038 


0.012 



TABLE II: Differential cross sections for the p+ He elastic 
scattering at 71 MeV/nucleon. The A# c .m. denotes the bin 
width. The A^fe denotes the quadratic sum of the statistical 
and systematic uncertainties. 



Oc.m. (de{ 


;) A8 c .m. (deg) 


Ay 


AAy 


37.1 


2.5 


0.242 


0.069 


44.6 


5.0 


0.021 


0.089 


54.6 


5.0 


-0.016 


0.135 


64.8 


5.0 


0.11 


0.18 


74.3 


5.0 


-0.27 


0.27 



TABLE III: Analyzing powers for the p+ He elastic scattering 
at 71 MeV/nucleon. AA y denotes the statistical uncertainty. 
Note that there is an additional scale error of 19% resulting 
from the uncertainty in the target polarization. The A^c.m. 
denotes the bin width. 
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FIG. 7: (Color online) Differential cross sections and analyz- 
ing powers of the p+ 4 He at 72 MeV (open squares: Ref. [23]), 
of the p+ 6 Li at 72 MeV (open triangles: Ref. [24]), and of 
the p+ 6 He at 71 MeV (open circles: Ref. [22], closed circles: 
present work). 



Ref. [22] in an overlapping angular region of 9 C m = 40 
50°. We extended the data to the backward angles of 
#c.m. ~ 90°. It is found that the da/dil of p+ 6 He are 
almost identical with those of p+ 6 Li at # c . m . = 20-90°, 
while they have a steeper angular dependence than those 
of p+ 4 Hc. In good contrast to the similarity found in 
da/dil, A y data are widely different between p+ 6 He and 
p+ 6 Li. The A y of p+ 6 Li increase as a function of the 
scattering angle in an angular region of 9 c . m , = 40-70° 
and take large positive values. This behavior is com- 
monly seen in proton elastic scattering from stable nu- 
clei at the present energy region [25]. Unlike this global 
trend, A y of p+ 6 Re decreases in 6> c . m . = 37-55°, which is 
rather similar to those of p+ 4 Hc. While the large error 
bars prevent us from observing the difference between A y 
of p+ 6 Hc and of p+ 4 Hc, it is clearly seen that the angu- 
lar distribution of the A y in p+ 6 Hc deviates from that of with 
p+ 6 Li. 



IV. PHENOMENOLOGICAL OPTICAL MODEL 
ANALYSIS 

A. Optical potential fitting 

The aim of this section is to extract the gross character- 
istics of the spin-orbit interaction between a proton and 



6 He. For this purpose, we determined the optical model 
potential that reproduces the experimental data of both 
differential cross sections and analyzing powers. The op- 
tical model potential obtained in this phcnomenological 
approach will be compared with the semi-microscopic cal- 
culations in Section V. 

We adopted a standard Woods-Saxon optical potential 
with a spin-orbit term of the Thomas form: 



C^om(-R) 



- V fr(R)-iW fi(R) 

+ Ha id W d ^-f id (R) 
dR 

+ V s ^^ R f s {R)L-a p + V c {R) (1) 



MR) 



1 + exp 



R - r^A 1 '* 



(2) 



r, i, id, or s). 



Here, R is the relative coordinate between a proton and 
a 6 He particle (see Fig. 12 (b)), L = Rx (-if£V R ) is 
the associated angular momentum, and <r p is the Pauli 
spin operator of the proton. The subscripts r, i, id, and 
s denote real, volume imaginary, surface imaginary, and 
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spin-orbit, respectively. Vc is the Coulomb potential of 
uniformly charged sphere with a radius of r^A 1 / 3 fm 
(r oc =1.3 fm). 

The search procedure for the best-fit potential param- 
eters was made in two steps: first the parameters of the 
central term were found by minimizing the x 2 values of 
da/dfl, and second the parameters of the spin-orbit term 
by fitting A y . These two steps were iterated alternately 
until convergence was achieved. Such a procedure is fea- 
sible since the contribution of the spin-orbit potential to 
da/dfl is much smaller than those of the central terms. In 
the fitting, we used the data in Ref. [22] and the present 
ones. Uncertainties of da/dfl smaller than 10% were ar- 
tificially set to 10% in order to avoid trapping in an un- 
physical local \ 2 minimum. The fitting was carried out 
using the ECIS79 code [26]. A set of parameters for the 
p+ 6 Li elastic scattering at 72 MeV/nuclcon [24], labeled 
as Set-A in Table IV, was used as the initial values in the 
search of the p- 6 He potential parameters. 

The parameters obtained for the p + 6 He elastic scat- 
tering are labeled as Set-B in Table IV. The reduced x 2 
values for da/dtt and A y were 0.95 and 0.96, respectively. 
Uncertainties of the parameters of the spin-orbit poten- 
tial, ros , a s j and V s , are evaluated in the following man- 
ner. Figure 8 shows the contour map of the deviation of 
X 2 value for A y from that calculated by the Set-B (as in- 
dicated by the point-P), &X A , on the two-dimensional 
plane of ro s and a s after projecting with optimized V s 
at each point of the plane. In the figure, a simultane- 
ous confidence region for ro s and a s is presented by the 
solid contour indicating Ax\ = 1. In this region, the 
optimum V s ranges between 1.15 MeV (at the point-Q) 
and 2.82 MeV (at the point-R). In the ro s -a s -V s space, 
a surface that = 1 touches planes that are ex- 

pressed by r 0s = 1.29 ± 0.13 fm, a s = 0.76 ± 0.17 fm, and 
V s = 2.02 ± 0.87 MeV, which gives a rough estimation 
for uncertainties of the parameters. 



B. Characteristics of spin-orbit potential 

In Fig. 9, the results of calculations of the observables 
made with the optical potentials of Set-A, -B, and -C in 
Table IV are shown together with the experimental data. 
Set-C was taken from Ref. [27], where a phenomenologi- 
cal optical model potential that reproduced only the pre- 
vious da/dVL data of the p+ 6 He at 71 MeV/nucleon [22] 
was reported. The radial dependences of the p- 6 He opti- 
cal potentials (Set-B and Set-C) are shown in Fig. 10 by 
solid and dashed lines, respectively. 

The calculation with the potential Set-C reasonably 
reproduces the present da/dQ data, whereas it largely 
deviates from the A y data at 9 c .m. ^ 40°. It should be 
noted that the A y data were unavailable when the poten- 
tial Set-C was sought. The calculation with the potential 
Set-B reproduces both da/dQ and A y over whole angular 
region except for the most backward data point of A y . 
Similarity of the da/dQ. calculated with Set-B and Set- 
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Os [ fm ] 

FIG. 8: Contour map of the Ax\ values (see the text for the 
definition) on the two-dimensional plane of ro a and a a . The 
solid, dashed, and dotted curves indicate Ax3i B =l, 3, and 
5, respectively. The point that gives the best-fit parameters, 
Set-B in Table IV, is indicated by the point-P. See the text 
for the points-Q and -R. 
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FIG. 9: (Color online) Differential cross sections and analyz- 
ing powers calculated by the phenomenological optical poten- 
tials are shown together with the experimental data. The dot- 
dashed curves denote calculations of the Set-A in Table IV, 
the solid curves those of set-B, and the dashed curves those 
of Set-C. Solid circles are present data and open circles are 
from Ref. [22]. 
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TABLE IV: Parameters of the optical potentials for p+ 6 Li at 72 MeV/nucleon [24] and p+ 6 He at 71 MeV/nucleon ([27] and 
present work). 
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Set-B 


p+ 6 He (Present) 
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0.76 


Set-C 


p+ G He [27] 


30.00 
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FIG. 10: (Color online) Radial dependences of the phe- 
nomenological optical potential (Set-B and Set-C in Ta- 
ble IV). Solid and dashed curves in the upper panel represent 
the real and imaginary parts of the central term. The lower 
panel shows the spin-orbit potential. 



C potentials originates from that of the central terms as 
shown in the upper panel of Fig. 10. The reliability of 
the potential obtained in the present work is supported 
by the fact that two independent analyses yielded the 
similar results for the central terms. In contrast to the 
central terms, the spin-orbit terms of these two potentials 
arc quite different, resulting in a large difference in A y as 
shown in Fig. 9. Note that the present data are sensitive 
to the optical potential in a region of R >1.5 fm. The 
spin-orbit potential of Set-B is much shallower than that 
of Set-C at R < 2.8 fm, while it is deeper at larger radii. 
This is due to the small value of V s and large values of 
ro s and a s of Set-B compared with those of Set-C. The 
phenomenological optical model analysis suggests that 
the A y data can be reproduced only with a shallow and 
long-ranged spin-orbit potential. 



The parameters of our spin-orbit potential are com- 
pared with those of neighboring even-even stable nuclei 
and with global potentials in Table V. Phenomenological 
optical potentials for the p+ 16 at 65 MeV and p+ 12 C 
at 16-40 MeV are taken from Ref. [25] and Refs. [28, 29], 
respectively. In addition to these local potentials, we 
also examined the parameters of global optical poten- 
tials: CH89 [4] and Koning-Delaroche (KD) [5], of which 
applicable ranges arc E = 10-65 MeV, A = 40-209 and 
E = 0.001-200 MeV, A = 24-209, respectively. While 
they are constructed for nuclei heavier than 6 He, it is 
worthwhile comparing them, since the mass-number de- 
pendence of the parameters is relatively small. For ex- 
ample, the mass-number dependence appears only in tq s 
in the case of CH89 [4] as: 

V s = 5.9(1) MeV, 
r 0s = 1.34(3) - 1.2(1)A~ 1/3 fm, 
a s = 0.63(2) fm. 

Table V includes the parameters of these potentials for 
the nuclei within the applicable range. Incident energies 
of E = 65 MeV and E =71 MeV were assumed here for 
CH89 and KD, respectively. 

Firstly, we focus on tq s and a s to discuss the radial 
shape of the spin-orbit potential. Combination of differ- 
ent values of tq s and a s can provide similar results of 
Ay since the observable is sensitive to the surface region 
of the spin-orbit potential. We thus compare these pa- 
rameters on the two-dimensional plane of ro s and a s as 
shown in Fig. 11. Parameters for the stable nuclei are 
mostly distributed in a region of ro s = 0.8-1.1 fm and 
a s = 0.5-0.6 fm, whereas that for 6 He is located in the 
upper right side of the figure. These large ro s and/or 
a s values indicate that the spin-orbit potential between 
a proton and a 6 He has a long-ranged nature compared 
with those for stable nuclei. The depth parameter V s 
was also compared with the global systematics. The V s 
value of p- 6 He potential was found to be 2.02 MeV for 
the best-fit potential (Set-B) and ranges between 1.15 
and 2.82 MeV in the simultaneous confidence region for 
ro s and a s . On the other hand, those of stable nuclei 
are mostly distributed around 5 MeV as shown in Ta- 
ble V. Comparing these values, the depth parameter of 
the spin-orbit potential between a proton and a 6 He is 
found to be much smaller than those of stable nuclei. 

The phenomenological analysis indicates that the spin- 
orbit potential between a proton and 6 Hc is characterized 
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V s I Xvlc V J 


ros (fm) 


ds (fm) 


p+ b Ho, E = 71 MeV (Set-B) 


2.02 


1.29 


0.76 


, 12ri p An MeV T9S1 

J>T Vy, Hi — 1V16 V [£o\ 


D.lo 


1 1 HQ 


U.01 ( 


p+ 12 C, E = 16-40 MeV [29] 


6.4 


1.00 


0.575 


p+ 16 0, £ = 65 MeV [25] 


5.793 


1.057 


0.5807 


CH89, E = 65 MeV, 4 =40-209 [4] 
KD, E = 71 MeV, 4 =24-209 [5] 


5.9±0.1 
4.369-4.822 


0.99-1.14 
0.961-1.076 


0.63±0.02 
0.59 



TABLE V: Parameters of the spin-orbit term of phenomenological and global optical potentials. 
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FIG. 11: Two-dimensional distribution of ro s and a a of 
phenomenological OM potentials for p+ He (closed circle), 
p+ 12 C (closed triangle: Ref. [29], closed diamond: [28]), and 
p+ 16 (closed square: [25]). The solid contour indicates the 
simultaneous confidence region for the ro a and a a values for 
the p+ 6 He as displayed in Fig. 8. Parameters of global OM 
potentials [4, 5] are also shown by solid lines which represent 
the ^-dependence of ro s . 



by large r$ s / a s and small V s values yielding shallow and 
long-ranged radial dependence. Intuitively, these char- 
acteristics can be understood from the diffused density 
distribution of 6 He. However, its microscopic origin can 
not be clarified by the phenomenological approach. To 
examine the microscopic origin of the characteristics of 
the p- e He interaction, microscopic and semi-microscopic 
analyses are required. Section V describes one of such 
analyses based on a cluster folding model for 6 He. 



V. SEMI-MICROSCOPIC ANALYSES 

In this section, we examine two kinds of the folding po- 
tential, the cluster folding (CF) and the nucleon folding 



(NF) ones. They are compared with the phenomeno- 
logical optical model (OM) potential determined in the 
preceding section. The results of calculations of observ- 
ables made by these potentials are compared with the 
experimental data. 

In the CF potential, we adopt the ann cluster model 
for 6 He and fold interactions between the proton and the 
valence neutrons, V pn , with the neutron density in 6 He 
and those between the proton and the a core, V pa , with 
the a density in 6 He. In the NF potential, we decompose 
the a core into two neutrons and two protons and fold 
the interactions between the incident proton and the four 
neutrons, V pn , with the neutron density in 6 He and those 
between the incident proton and two target protons, V pp , 
with the proton density in 6 He. 

The detailed expressions of such folding potentials are 
given in the following subsection, where the Coulomb in- 
teraction is considered in the p-p and p-a interactions 
respectively when compared with the corresponding scat- 
tering data but finally it is assumed to act between the 
proton and the 6 He target with roc = 1.400 fm [27]. 

A. Folding potentials 

Denoting two valence neutrons by ri\ and n 2 , the CF 
potential C/cf is given as 

U C f = J I/ p „ 1 p^ F (r 1 )dr 1 + J V pn2 p° F {r 2 )dr 2 

+ J V paP ^ F (r a )dr a , (3) 

where ri, r 2 , and r a are the position vectors of n±, n 2 , 
and the a core from the center of mass of 6 He, respec- 
tively. The neutron and a densities, p„ and p„ F , are 
calculated by the ann cluster model for 6 He [30, 31], 
where the condition ri + r 2 + 4r a = is considered as 
usual. 

In the present work, we specify the potentials in the 
right hand side of Eq. (3) by the central plus spin-orbit 
(LS) type: 

Vpn, = C(l r P». I) + Vpndrpn, ' K + ^J, 

where i = 1 , 2 and 

V pa = V p ° a (\r pa \) + V p L *(\r pa \)£ pa ■ a p . (4) 
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FIG. 12: (Color online) Coordinate systems for the cluster 
folding model. 



Here, r pnil r pn2 , and r pa are denned in Fig. 12 (a), and 

^pni Tpni X ( ith\/ pni ) s etc. 

In the following, we transform the set of coordinates 
(fprai j r pn 2 ' r pa) to that of (£, £, Jf2), which are defined in 
Fig. 12 (b), to describe the angular momenta l pni and 
t pa in terms of L. The transformation is 



pet 



= R 
= -B + 



1 



(5) 



and consequently 









-v £ , 






6 








^ pa 


2^ 






(G) 



These relations lead to, for example, 



V C )] • (7) 



Here, and can be neglected, because these are 
the momenta for the internal degrees of freedom of 6 He 
and their expectation values are zero for a spherically 
symmetric nucleus [32]. Using |£ + i£ = — r x , we get 



£ pni =j[L-nx (-iftVfl)] , (8) 



which is independent of the special choice of the 6 He 
internal coordinates, £ and £. To i, ri can contribute 
by its component along the R direction [32] , then 



l -L(l nR ) 



(9) 



Similar expressions are obtained for £ pn2 and £ pa . Setting 
5(°Vu + <T n 2 ) = and considering other quantities to 
appear in symmetric manners on 1 and 2, we obtain the 
p- 6 He potential as 



U CF = C/ (R) + Uti (R) L-a v , 



(10) 



with 



U£ F (R) = 2 |^(|n-Jl|)pS P (ri)dri 
+ y F p a (|r a -H|)p2 F Mdr Q 



(11) 



and 



rLS " r a R 



V£(\r a -R\) 1 



i? 2 



In a way similar to the above development, we get the 
NF model potential C/nf- In this case, the relative coor- 
dinates between the incident proton and six nucleons in 
the 6 He nucleus are transformed to the proton- 6 He rela- 
tive coordinate R and a set of five independent internal 
coordinates of 6 He. The obtained J/nf, which is indepen- 
dent on the choice of the set of the internal coordinates, 
is written as 



with 



C/nf - U$ F (R) + U^i(R) L-rr p , (13) 

U* F (R) = 2 J l^drx-flDpf (n)dn 

+ 4 Jv p ° n (\r 2 ~R\)p* F (r 2 )dr 2 (14) 



and 

U?i(R) 



i [vh S (\r 1 -R\){l-^) P "/(r 1 )dr 1 



pp v 

rLS/ 



1 - 



R 2 
r 2 R 



.NF 



where /?^ F and p p F denote point neutron and proton den- 
sities, respectively. 



B. Numerical evaluation of p- He potentials 

To evaluate the p- 6 He folding potentials as specified 
in the preceding section, we have to fix the following 
elements; the p-a interaction V pai the p-p and p-n in- 
teractions V pp and V pn , and the densities in 6 He, p a , p P 
and p n . These are discussed in the following subsections, 
respectively. 



p-a interactions 



For V pa 



used in the CF potential, we assume the stan- 
dard WS potential such as given in Eq. (1). The pa- 
rameters involved are searched so as to fit the data of 
dcr/dVt and A y in the p+a scattering at 72 McV/nucleon 
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[23]. Particular attention was given to reproducing the 
observables in the forward angular region, since overall 
agreements with the data are not found in spite of the 
careful search of the parameters. Two typical parameter 
sets, with and without the volume absorption term, are 
labeled as Set-1 and Set- 2 in Table VI. The results of cal- 
culations made with these potentials arc compared with 
the data in Fig. 13, where the solid and dashed lines show 
those by Set-1 and Set-2 potentials, respectively. Both 
calculations describe the data up to 9 w 100° but do 
not reproduce those at backward angles, 9 > 120°. Such 
discrepancies between the calculated results and the mea- 
sured data at the backward angles suggest participation 
of contributions of other reaction mechanisms, such as 
knock-on type exchange scattering of the proton with 
target nucleons. Such possible extra mechanisms will 
be disregarded at present since we are concerned with 
the p-a one-body potential. In our CF calculations, we 
adopt the potentials with the above parameter sets as 
V pa . However, the validity of the CF potential thus ob- 
tained is limited to forward scattering angles, a low mo- 
mentum transfer region, of p+ e He scattering. The real 
and imaginary parts of V pa and the real part of V^f for 
the above parameter sets are displayed in the upper and 
lower panels of Fig. 14. Although Set-1 (dashed) and 
Set-2 (solid) potentials have rather different r pa depen- 
dence, as shown later, this difference is moderated in the 
folding procedure so yielding similar CF potentials. 



2. p-p and p-n interactions 
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FIG. 13: (Color online) Angular distribution of the 
cross section and A y for the p+ 4 He elastic scattering at 
72 MeV/nucleon. The solid and dashed lines are the optical 
model calculations with Set-1 and Set-2 parameters, respec- 
tively. The experimental data are taken from Ref. [23]. 



For V pn and V pp used in the CF and NF potentials, 
we adopt the complex effective interaction, CEG [33-35] , 
where the nuclear force [36] is modified by the medium 
effect which takes account of the virtual excitation of nu- 
cleons of the nuclear matter up to 10kp by the ^-matrix 
theory. The nuclear force is composed of Gaussian form 
factors and the parameters contained are adjusted to sim- 
ulate the matrix elements of the Hamada- Johnston po- 
tential [37]. The CEG interaction has been successful in 
reproducing dcr/dil and A y measured for the proton elas- 
tic scattering by many nuclei in a wide incident energy 
range, E p — 20-200 MeV, in the framework of the fold- 
ing model [33-35] . It has been shown that the imaginary 
part of the folding potential given by the CEG interac- 
tion is slightly too large to reproduce experimental N-A 
scattering [33, 35]. In the present calculation, therefore, 
we adopt the normalizing factor Ni = 0.7 for the imagi- 
nary part of the CEG interaction. However, calculations 
with Nj = 1.0 do not give an essential change to the 
results. 



3. Densities of a, p and n in 6 He 

The densities p^ F and p^ F for the CF calcula- 
tion arc obtained by applying the Gaussian expansion 



method [30, 31] to the ann cluster model of 6 He. This 
method has been successful in describing structures of 
various few-body systems as well as 6 He [30, 31]. As 
for the n-n interaction, we choose AV8' interaction [38]. 
It is reasonable to use a bare (free space) n-n interac- 
tion between the two valence neutrons in 6 He as they 
are dominantly in a region of low density. As for the 
a-n interaction, we employ the effective a-n potential in 
Ref. [39], which was designed to reproduce well the low- 
lying states and low-energy-scattering phase shifts of the 
a-n system. The depth of the a-n potential is modi- 
fied slightly to adjust the ground-state binding energy of 
6 He to the empirical value. In Fig. 15(a), the densities 
obtained are shown as functions of r, the distance from 
the center of mass of 6 He, where is localized in a 
relatively narrow region around the center, while p^ F is 
spread widely. 

The NF calculation depends on the assumptions made 
for the densities of the two protons and four neutrons 
in 6 He as well as those made for the p-p and p-n in- 
teractions [18]. At present, to see the essential role of 
clustering the four nucleons into the a-particle core, we 
use the densities of the proton and the neutron in the 
a obtained by decomposing the density of the point a, 
Pq F , to the densities of the constituent nucleons with a 
one-range Gaussian form factor with range 1.40 fm. The 
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TABLE VI: Parameters for the optical potentials for p+ 4 He at 72 MeV/nucleon. 
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FIG. 14: (Color online) Optical potentials for the p+ 4 He elas- 
tic scattering at 72 MeV/nucleon. The thick dashed (solid) 
lines are for the real part of Set-1 (Set-2) potential and the 
thin dashed (solid) lines are for the imaginary part. 



total nucleon densities of 6 He, p p F , and p^ F , where the 
latter includes the contribution of the valence neutrons, 
are displayed in Fig. 15(b). The neutron density p^ F has 
a longer tail than the proton one p p F due to the presence 
of the valence neutrons. In Refs. [40, 41] the nucleon den- 
sities of 6 He were calculated in a more sophisticated way. 
They produced densities similar to the present ones for 
the protons and neutrons. These two kinds of nucleon 
densities provide similar results in the NF calculation of 
dcr/dVt and A y of the p+ 6 He scattering. Thus, in the fol- 
lowing, wc will discuss J7nf as formed using the densities 
shown in Fig. 15(b). 



4- p- 6 He folding potentials 

In Fig. 16, the resultant p- 6 He potentials, Ucf and 
£7nF; arc compared with each other as well as with the 
optical model potential Uom- The CF potentials cal- 
culated by the two sets of V pa in Tabic VI, say Ucf-i 
and Ucf-2, are shown by long-dashed and solid lines in 
Fig. 16(a), respectively. The folding procedure gives sim- 
ilar results in both cases. The contribution of V pn is dis- 
played by short-dashed lines in the figure, which is found 
to be mostly small. Especially, in the spin-orbit poten- 
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FIG. 15: (Color online) The densities of a, p, and n in 6 He 
used in folding models. The p„ F and p^ F in panel (a) are 
used in the CF calculation. The p^ F and p^ F in panel (b) 
are used in the NF calculation. All densities are normalized 
as 47r f p x r 2 dr — 1, where x = a, n, and p. 



tial, the contribution from V pn is one order of magnitude 
smaller than that from V pa . The main contribution to 
Uqf arises from the interaction V pa except for the cen- 
tral real potential at R >3 fm, which is dominated by 
the Vpn contribution. This is supposed to be the reflec- 
tion of the extended neutron density shown in Fig. 15 
and produce significant contributions to the observables 
as discussed later. 

In Fig. 16(b), Ucf due to Set-2 of V pa , J7nf, and Uom 
are shown by solid, dot-dashed, and short-dashed lines, 
respectively. First we consider the central part of the 
potentials. For small R, the real part of Uq F is deeper 
than those of C/ NF and C/ OM , while for large R, C/ CF is 
shallower than other two. In the imaginary part, the 
magnitude of Uq is much bigger than those of the other 
two potentials. This will compensate the deficiency of the 
real part of Uq F at large R, for example in the calculation 
of the cross section. On the other hand, the spin-orbit 
part of Uom has larger magnitude for R > 2 fm and 
thus has a longer range compared with those of other 
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FIG. 16: (Color online) Potentials between proton and 6 He. (a): The long-dashed (solid) lines are the CF calculation with Set-1 
(Set-2) parameters for V pa - The short-dashed lines show the contribution of the valence two neutrons to the CF potential, (b): 
The solid lines are the CF calculation with Set-2 parameters, the dot-dashed lines are the NF one and the short-dashed lines 
are the phenomenological optical potential. The CF spin-orbit potential corrected by AC/ls term is shown by the long-dashed 
line (see text for detail). 



two potentials. Such long-range nature of the spin-orbit 
interaction is a characteristic feature of the spin-orbit 
part of Uom as described in Section IV. This is discussed 
later in more detail with relation to A v . 



C. Comparison between experiments and 
calculations in p+ 6 He scattering 

In the following, the da/dQ, and A y for p- 6 He elas- 
tic scattering calculated using Uqf, E^nf> and Uom are 
compared with the data taken at 71 MeV/nucleon. In 
Fig. 17(a) the results obtained using the two CF poten- 
tials, C/cF-i and Ucf-2, are shown by long-dashed and 
solid lines, respectively. Both results are very similar 
to each other and well describe the data of da/dtt, ex- 
cept for large angles where the calculations overestimates 
the data by small amounts. The calculations also de- 
scribe the angular dependence of the measured A y up 
to 9 ~ 55°. These successes basically support the CF 
potential as a reasonable description of the scattering. 
The discrepancies at large angles, i.e. a large momentum 
transfer region, may be related to the limitation of the 
validity of V pa used in the folding, as discussed in the 
subsection VB. 

In Fig. 17(a), the results of the calculation made using 
the NF potential are shown by dot-dashed lines. They 
do not reproduce the data well. The calculation gives an 
deep valley around 8 ~ 53° in the angular distribution 
of da/dil and a large positive peak at the correspond- 
ing angle of the A y angular distribution. These features 



do not exist in the data. Since the present nucleon den- 
sities originated from the CF model ones, the essential 
difference between the CF and NF potentials will be pro- 
duced by the use of the different interactions. Thus, the 
CF calculation will owe its successes to the inclusion of 
the characteristics of the realistic p-a interaction into the 
p- 6 He potential. 

It is interesting to examine if the a core in 6 He is some- 
what diffused compared with a free a-particle, due to the 
interactions from the valence neutrons. For that purpose, 
we increased the radius and diffuseness parameters, r 
and a, in V pa potential as r$ to 1.1 rg and a to a + 0.1 
fm. The depth parameters were changed to keep constant 
the values of the corresponding volume integrals. The ef- 
fect of this change is shown by the short-dashed lines in 
Fig. 17(a), where reproduction of the data is improved 
somewhat, especially in da/dfl at large angles. 

In Fig. 17(b), the contributions of the valence neu- 
trons are demonstrated for the CF-2 calculation. As is 
speculated from the analyses of the form factors of the 
potential in Fig. 16(a), the dominant contribution to the 
observables in the CF calculation arises from V pa dis- 
played by the dot-dashed lines in Fig. 17(b). However, 
the valence neutrons produce indispensable corrections 
to the observables. That is, the pn central interaction 
decreases da/dQ at large angles, giving remarkable im- 
provements of the agreement with the data as shown by 
the dashed lines. The pn interaction also contributes to 
A y by a considerable amount through the central part. 
A detailed examination of the calculation revealed that 
such corrections were due to the V pn part of the folding 
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FIG. 17: (Color online) Angular distribution of the cross section and A v for the p+ 6 He elastic scattering at 71 MeV/nucleon. 
The experimental data are denoted by circles (present) and squares (Ref. [22]). (a): The long-dashed (solid) lines are the CF 
calculation in which Set-1 (Set-2) parameters are used for V pa - The dot-dashed lines are the NF calculation, (b): The dashed 
lines and the solid ones include the V p „ interaction, where the formers neglect the spin-orbit part of V pn - The dash-dotted lines 
include only V pa interaction, (c): The solid lines are the CF calculation with Set-2 parameters for V pa and the short-dashed 
lines are the OM calculation. The long-dashed lines are the CF calculation with Set-2 parameters for V pa where A£/ls is added 
(see text for detail). 
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central potential in a R region between R =2 fm and 4 fm 
(see Fig. 16(a)). The spin-orbit part of V pn gives almost 
no effect to the observables as shown by the solid lines in 
Fig. 17(b). This is consistent with the result shown by 
Crespo et al. [42] in a study at higher incident energy. 

In Fig. 17(c), we compared the results of the CF-2 cal- 
culation (solid lines) with those of the OM calculation 
in the preceding section (long-dashed lines) as well as 
with the data. In the optical model analysis, the experi- 
mental data can be reproduced only with a shallow and 
long-ranged spin-orbit potential. Compared with this po- 
tential, the spin-orbit part of CF-2 potential has a shorter 
range as displayed by a solid line in the lower panel of 
Fig. 16(b). To investigate the role of the long tail in 
the spin-orbit potential, we calculated the observables 
by adding a weak but long-range spin-orbit interaction 
AC/ls to the CF interaction. This correction is assumed 
to be the Thomas type as 



AC/ L s OR) 



2 d 



VaM R dR 



1 + cxp 



{(R 



6 1/3 r add )/a add 

(16) 

For simplicity, we adopt v add = 1 MeV, r add = 1.5 fm, 
and a add = 0.7 fm where the large magnitudes of r add 
and a add are consistent with the characteristics of the 
magnitudes of ro s and a s of the OM potential discussed 
in Sec. IV. The calculated observables are displayed in 
Fig. 17(c) by short-dashed lines, where da/dVl is little 
affected but A y receives a drastic change, i.e. the angu- 
lar distribution of A y is now similar to that by the OM 
calculation in a global sense showing qualitative improve- 
ments in comparison with the data. To see the contri- 
bution of AC/ls to the potential, we plot + AC/ls 
in Fig. 16(b) by long-dashed lines, where the new spin- 
orbit potential becomes very close to that of the OM 
potential at R > 2.5 fm. It is indicated that the long 
tail of the spin-orbit potential is particularly important 
in reproducing the angular distribution of A y , while its 
microscopic origin is still to be investigated. When some 
corrections which increase the range of the spin-orbit in- 
teraction are found, they will be effective for improving 
the CF calculation. 



VI. MICROSCOPIC MODEL ANALYSES 

In this section, we will describe the theoretical analysis 
of the present data by a microscopic model developed in 
Ref. [7]. In this model, one can predict the scattering 
observables such as cross sections and analyzing pow- 
ers with one run of the relevant code (DWBA98) with 
no adjustable parameter. Complete details as well as 
many examples of use of this coordinate space micro- 
scopic model approach are to be found in the review [7] . 
Use of the complex, non-local, nucleon-nucleus optical 
potentials defined in that way, without localization of 
the exchange amplitudes, has given predictions of differ- 
ential cross sections and spin observables that are in good 
agreement with data from many nuclei ( 3 Hc to 238 U) 



and for a wide range of energies (40 to 300 MeV) . Cru- 
cial to that success is the use of effective nucleon-nucleon 
(NN) interactions built upon NN g-matrices. The ef- 
fective NN interactions are complex, energy and density 
dependent, admixtures of Yukawa functions. They have 
central, two-nucleon tensor and two-nucleon spin-orbit 
character. The NA optical potentials result from fold- 
ing those effective interactions with the one-body density 
matrix elements (OBDME) of the ground state in the tar- 
get nucleus. Antisymmetrization of the projectile with 
all target nucleons leads to exchange amplitudes, making 
the microscopic optical potential non-local. For brevity, 
the optical potentials that result are called g-folding po- 
tentials. Another application has been in the prediction 
of integral observables of elastic scattering of both pro- 
tons and neutrons, with equal success [43]. Thus, the 
method is known now to give good predictions of both 
angular-dependent and integral observables. 
_j It is important to note that the level of agreement with 
.data in the ^-folding approach depends on the quality of 
the structure model that is used. Due to the character 
of the hadron force, proton scattering is preferentially 
sensitive to the neutron matter distributions of nuclei; 
a sensitivity seen in a recent assessment, using proton 
clastic scattering, of diverse Skyrme-Hartree-Fock model 
structures for 208 Pb [44]. 



A. Structure of He used 

6 He is a two-neutron halo nucleus and has been de- 
scribed well by shell model calculations. In calculation of 
the (/-folding potential for protons interacting with 6 He, 
a complete (0+2+4)fiu; shell model calculation has been 
made to specify the ground state OBDME. Essentially 
they are the occupation numbers which define the mat- 
ter densities of the nucleus. 

In the present study, we assume three sets of the single- 
nucleon (SN) wave functions for 6 He. One is the oscil- 
lator wave functions with an oscillator length of 2.0 fm 
(HO set). However, a neutron- halo character of 6 He can 
not be given by the oscillator wave function whatever os- 
cillator length is used as shown by the dashed curve in 
the lower panel of Fig. 18. Thus, we assume two sets of 
SN wave functions defined in Woods-Saxon (WS) poten- 
tials. One of them is obtained by taking the geometry 
of the potential from that found appropriate in Ref. [7], 
where electron form factors and proton scattering from 
6,7 Li are studied. That study provided a set of SN wave 
functions that we specify as WS nonhalo set since the p- 
shell nucleons were all reasonably bound. The extended 
neutron matter character of 6 Hc is found by choosing 
the binding energy of the halo-neutron orbits to give the 
single-neutron separation energy (1.8 MeV) to the lowest 
energy resonance in 5 He. The set of SN wave functions 
that result are specified as WS halo set. The associated 
density profile has the extensive neutron density coming 
from the halo. Density profiles given by the various sets 
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of SN wave functions are shown in Fig. 18. The dashed, 
solid, and dot-dashed curves show the density distribu- 
tions of HO, WS halo, and WS nonhalo sets, respectively. 
The difference between proton distributions of WS halo 
and WS nonhalo sets can not be seen. 
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FIG. 18: The (model) proton and neutron densities are shown 
in the upper and lower panels, respectively. The solid, dot- 
dashed, and dashed curves represent those obtained with WS 
halo, WS nonhalo, and HO sets, repectively. 



of SN wave functions. This calculation does not give a 
satisfactory result; especially in the case of the analyzing 
power. The solid curves show results found using WS 
halo set while those depicted by the dot-dashed curves 
are those found with WS nonhalo set. Of these the halo 
description gives the better match to data especially at 
the larger scattering angles. This result is consistent with 
the findings from analyses of lower energy scattering data 
at 40.9 MeV/nucleon [45] and at 24.5 MeV/nucleon [46]. 
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Use of WS halo set in analyses of 40.9 MeV/nucleon 
data [45] gave a value of 406 mb for the reaction cross 
section, which is in good agreement with the measured 
value. Additional evidence for WS halo set is given by 
the root mean square (r.m.s.) radius of the matter dis- 
tribution, which is most sensitive to characteristics of 
the outer surface of a nucleus. Using WS nonhalo set of 
the SN wave functions gave an r.m.s. radius for 6 He of 
2.30 fm, which is much smaller than the expected value 
of 2.54 fm. On the other hand, using WS halo set gave 
an r.m.s. radius for 6 He of 2.59 fm in good agreement 
with that expectation. 



B. Differential cross sections and analyzing powers 

The cross sections and analyzing powers for the p+ 6 Hc 
elastic scattering at 71 MeV/nucleon are shown in the 
top and bottom panels of Fig. 19, respectively. The cal- 
culated results shown therein by the dashed lines were 
found using the ^-folding potential obtained with HO set 



FIG. 19: Differential cross sections and analyzing powers of 
the p+ 6 He elastic scattering at 71 MeV/nucleon (open cir- 
cles: Ref. [22], closed circles: present work). Three curves 
are results of g-matrix folding calculation with 6 He densities 
presented in Fig. 18. 



In fact, the differential cross sections calculated with 
WS halo set match the data so well that one does not 
need to contemplate any adjustment. However, the story 
is not so simple when one also considers the analyzing 
power data. At forward scattering angles, both WS sets 
reasonably match the data. But neither WS result pro- 
duces the distinctive trend of small values found at larger 
scattering angles. Nonetheless the best result is that 
found on using WS halo set of SN functions. Given that 
the cross section values in the region of 60° to 90° is of an 
order of 0.1 mb/sr, the limitations in the present micro- 
scopic model formulation of the reaction dynamics may 
be the problem. 
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VII. SUMMARY 

The vector analyzing power has been measured for 
the elastic scattering of 6 He from polarized protons at 
71 MeV/nucleon to investigate the characteristics of the 
spin-orbit potential between the proton and the 6 He nu- 
cleus. Measurement of the polarization observable was 
realized in the Rl-beam experiment by using the newly 
constructed solid polarized proton target, which can be 
operated in a low magnetic field of 0.1 T and at high tem- 
perature of 100 K. The measured da/dCl of the p+ 6 He 
elastic scattering were almost identical to those of the 
p+ 6 Li. On the other hand, the A y were found to be 
largely different from those of the p+ 6 Li and rather sim- 
ilar to those of the p+ 4 He elastic scattering. 

To extract the gross feature of the spin-orbit interac- 
tion between a proton and 6 He, an optical model po- 
tential was determined phcnomcnologically by fitting the 
experimental data of da/dfl and A y . Compared with the 
global systematics of the potentials for stable nuclei, it is 
indicated that the spin-orbit potential for 6 He is charac- 
terized by a small value of V s and large values of ro s and 
a Sl namely by a shallow and long-ranged radial shape. 
Such characteristics might be the reflection of the dif- 
fused density of the neutron-rich 6 He nucleus. 

The cluster folding calculation was carried out to get 
a deeper insight into the optical potential, assuming the 
a+n+n cluster structure for 6 He. In addition, nucleon 
folding calculations were also performed by decomposing 
the a core into four nucleons. The experimental data 
could not be reproduced by the nucleon folding calcula- 
tion, whereas the arm cluster folding calculation gives 
the reasonable agreements with the data. Thus, this in- 
dicates that it is important to take into account of the 
a-clustcrization in the description of p+ 6 He elastic scat- 
tering. The cluster folding calculation shows that the 
dominant contribution to the p- 6 Hc potential arises from 
the interaction between the proton and the a core. Es- 
pecially, in the spin-orbit potential, the contribution of 
the interaction between the proton and valence neutrons 
was found to be much smaller than the a core contri- 
bution. However, the measured cross section at large 
angles can not be understood without the contribution 
from the scattering by the valence neutrons. Comparison 
of the phenomenological optical potential and the cluster 
folding one indicates that the long-range nature of the 
spin-orbit potential is important in reproducing the A y 
data at large angles. The microscopic origin of such a 
long tail is still to be investigated. 

The data were also compared with the predictions ob- 
tained from a fully microscopic g-folding model. Three 
sets of single nucleon wave functions were tried since 
other details of the calculation were predetermined. The 
model, which has been successful in analyzing p+ 6 He 
scattering cross sections in the past [44] , again gives good 
reproduction of the data in the present case when the 
bound state wave functions specify that 6 He has a neu- 
tron halo. However, the match to the data, in particular 



the analyzing power, is not perfect. This may indicate 
limitation of the structure model used and/or of unac- 
counted reaction mechanisms that influence the larger 
momentum transfer results. 

This work has demonstrated the capability of the solid 
polarized proton target in low magnetic field and high 
temperature to probe the new aspects of the reaction 
involving unstable nuclei. Future polarization studies of 
such kinds will provide us with valuable information on 
the reaction and structure of unstable nuclei. 
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Appendix A: Absolute measurement of target 
polarization 

In the case of conventional solid polarized targets, the 
NMR signal usually is related to the absolute magni- 
tude of the polarization by measuring the target po- 
larization under the state of thermal equilibrium (TE). 
However, measurement of the TE polarization is quite 
difficult in our target. The first reason for this is that 
the TE polarization is very small in a low magnetic 
field and at high temperature, since it is represented by 

-Pte = tanh (jfijp^j , where fi, B, and T are the magnetic 
moment of proton, the field strength, and the tempera- 
ture, respectively. The second reason is that the sensi- 
tivity of the present NMR system is not sufficiently high 
since the target design is optimized for scattering exper- 
iments. 

One of the simple methods to measure the absolute tar- 
get polarization would be the measurement of the spin- 
dependent asymmetry e = P y A y for the proton elastic 
scattering whose analyzing power A y are known. In the 
present study, we measured the spin-asymmetry for the 
p+ 4 He elastic scattering at 80 MeV/nucleon. The A y 
have already been measured by Togawa et al [47]. The 
use of the p+ 4 He scattering is profitable since we can 
measure the e with the same experimental setup as that 
for the p+ 6 He measurement only by changing settings 
of the fragment separator RIPS to produce a secondary 
4 He beam. The profile of the 4 He beam on the target 
was tuned to be almost same as that of the 6 He beam. 

Figure 20 shows A y of the p+ 4 He elastic scattering 
at 80 MeV/nucleon. The open circles represent the pre- 
vious data [47], while the closed ones show the present 
data whose magnitudes are scaled to the previous ones. 
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From the scaling factor, the average polarization dur- 
ing the p+ He measurement was determined to be P y = 
12.3 ± 2.4%. The relative uncertainty of the polarization 
APy/Py, which was 19% in the present work, resulted 
from the statistics of the p+ 4 He scattering events. Fu- 
ture development of the NMR system would be required 
for determining the absolute polarization more precisely 
without losing beam time. 
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FIG. 20: Analyzing powers of the p+ 4 He elastic scattering 
at 80 MeV/nucleon. Closed circles indicate the present data 
where P y = 12.3% is assumed. Open circles represent the 
reference data taken from Ref. [471 . 
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